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Abstract 

In a semiclassical context we investigate the Zitterbewegung of relativistic par- 
ticles with spin 1/2 moving in external fields. It is shown that the analogue of 
Zitterbewegung for general observables can be removed to arbitrary order in h by 
projecting to dynamically almost invariant subspaces of the quantum mechanical 
Hilbert space which are associated with particles and anti-particles. This not only 
allows to identify observables with a semiclassical meaning, but also to recover com- 
bined classical dynamics for the translational and spin degrees of freedom. Finally, 
we discuss properties of eigenspinors of a Dirac-Hamiltonian when these are projected 
to the almost invariant subspaces, including the phenomenon of quantum ergodicity. 
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1 Introduction 



When Dirac had introduced his relativistic quantum theory for spin-1/2 particles it be- 
came soon clear that despite its overwhelming success in producing the correct hydrogen 
spectrum, including fine structure, this theory was plagued with a number of apparent in- 
consistencies as, e.g., Klein's paradox |Kle29j . Schrodinger observed |Sch30j that the (free) 
time evolution of the naive position operator, that he had taken over from non-relativistic 
quantum mechanics, contained a contribution without any classical interpretation. Since 
this term was rapidly oscillating he introduced the notion of Zitterbewegung (trembling mo- 
tion). Schrodinger wanted to explain this effect as being caused by spin. Later, however, 
it became clear that all paradoxa stemmed from the coexistence of particles and anti- 
particles in Dirac's theory, and that only a quantum field theoretic description, predicting 
pair creation and annihilation, could cure these problems in a fully satisfactory manner. 

It nevertheless turned out that relativistic quantum mechanics was capable of describ- 
ing many important physical phenomena; e.g., in atomic physics it serves to explain spectra 
of large atoms and in nuclear physics it allows to construct models of nuclei. Within the 
context of this theory one can remove the Zitterbewegung of free particles by introducing 
modified position operators that are associated with either only particles or only anti- 
particles. These operators are free from the particle/anti-particle interferences that cause 
the Zitterbewegung; for details see e.g. |Tha92| . In the case of interactions (with exter- 
nal potentials) the Zitterbewegung cannot be exactly eliminated. It is, however, possible 
to devise an asymptotic construction, e.g. based on the Foldy-Wouthuysen transforma- 
tion. This involves a non-relativistic expansion, and when applied to the Dirac equation 
reproduces in leading order the Pauli equation. The genuinely relativistic coexistence of 
particles and anti-particles is therefore removed in a natural way. 

Similar constructions are possible in a semiclassical context while maintaining the rel- 
ativistic level of description. The decoupling of particles from anti-particles then proceeds 
asymptotically order by order in powers of the semiclassical parameter h. In this direc- 
tion several approaches have been developed recently, aiming at semiclassical expansions 
for scattering phases [BN99[ BR99 or at recovering the Thomas precession of the spin 
|SpoOO , see also |CorOH ITeu03j. In other semiclassical approaches trace formulae of the 



Gutzwiller type have been derived |BK98l IBK99j , see also |Kep03| . Here our goal is to 
introduce a general semiclassical framework for the Dirac equation, and in particular to 
identify a sufficiently large algebra of quantum mechanical observables that can be as- 
signed a classical analogue in a way that is consistent with the dynamics. And in contrast 
to previous studies we do not construct semiclassical (unitary) Foldy-Wouthuysen opera- 
tors but use the associated projection operators, thus avoiding ambiguities in the choice of 
the unitaries. 

The outline of this paper is as follows. In section [2] we recall the phenomenon of 
Zitterbewegung of free relativistic particles, emphasising that this effect can be removed 
through an introduction of suitably modified position operators. The following section is 
devoted to the construction of projection operators that semiclassically separate particles 
from anti-particles in the presence of an interaction with external potentials. Quantum 
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observables with a semiclassical meaning are introduced in section HJ and a characterisa- 
tion of such observables that maintain this property under their time evolution is provided. 
Section El contains a discussion of the classical dynamics emerging from the quantum dy- 
namics generated by a Dirac-Hamiltonian in the semiclassical limit. It is shown that a 
combined time evolution of the translational and the spin degrees of freedom arises in the 
form of a so-called skew-product. Finally, in section we discuss to what extent eigen- 
spinors of a Dirac-Hamiltonian can be associated with either particles or anti-particles. We 
also show that an ergodic behaviour of the combined classical dynamics implies quantum 
ergodicity, i.e. a semiclassical equidistribution on particle and anti-particle energy shells 
of eigenspinor-projections to the respective energy shells. 



2 Zitterbewegung for the free motion 

Since its early discovery by Schrodinger Zitterbewegung has mostly been discussed in the 
context of the Dirac equation for a free relativistic particle of mass m and spin 1/2, whose 
Hamiltonian reads 

H := -ikca ■ V + f3mc 2 . (2.1) 

The relations a^ai + otiotk = 25m and ak/3 + (3ak — (k,l = 1,2,3) defining the Dirac- 
algebra are realised by the 4x4 matrices 

oc= {l o ) and o 2 -1 2 ) ' 

where cr = (o"i, a 2 , cr 3 ) is the vector of Pauli spin matrices and 1„ denotes the n x n unit 
matrix. The operator (j2.1j) is essentially self-adjoint on the dense domain C^°(1R 3 ) <S> C 4 in 
the Hilbert space := L 2 (1R 3 ) £g> C 4 ; its spectrum is absolutely continuous and comprises 
of (— oo, — mc 2 ) U (mc 2 , oo), see e.g. [Tha92j. The time evolution (free motion) 

U (t) := e"^ ' 

hence is unitary. 

In his seminal paper Sch30j Schrodinger solved the free time evolution of the standard 
position operator x, whose components are defined as multiplication operators on a suitable 
domain in L 2 (IR 3 )®C 4 . His intention was to resolve the apparent paradox that the spectrum 
of the velocity operator resulting from the Heisenberg equation of motion for the standard 
position operator, x(t) = ca(t), consists of ±c, whereas the classical velocity of a free 
relativistic particle reads c 2 p/ \Jp 2 c 2 + m 2 c 4 and thus is smaller than c in magnitude. 
Integrating the equations of motion Schrodinger found 

x(t) = U (t)*xU (t) 

= x(0) + c 2 pH H + ^H \^ U ' 2] 
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The first two terms of this result exactly correspond to the respective classical dynamics 
of a free relativistic particle. The third term, however, contains the operator 



F := col - c 2 pH 1 , 

which is well defined since zero is not in the spectrum of H and thus Hq 1 is a bounded 
operator. The quantity F expresses the difference between the velocity operator and the 
quantisation of the classical velocity. This additional term introduces a rapidly oscillating 
time dependence and hence was named Zitterbewegung (trembling motion) by Schrodinger. 

After Schrodinger's work Sch30] the origin of the Zitterbewegung was traced back to 
the coexistence of particles and anti-particles in relativistic quantum mechanics. Therefore, 
this effect can be removed by projecting the standard position operator to the particle and 
anti-particle subspaces of J$?, respectively. To this end one introduces the projection 
operators 

P±:=Ul je ±\H \- 1 H ) , (2.3) 
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with \H \ := (#*F ) 1/2 , fulfilling P+P = and P+ + P = 1#. The time evolution of 
the projected position operators 

xt:=P*xF$ (2.4) 

is then given by (see |Tha92j ) 

x±(t) = x±(0) + c 2 P H^tP ± . 

It hence exactly corresponds to the respective classical expression; and this is true for both 
particles and anti-particles. The interpretation of this observation is obvious: Particles 
and anti-particles show noticeably different time evolutions and the Zitterbewegung in 
(j2.2j) represents the interference term [T ha92j 

ifj . 

p+^,p- — I — p-Z,p+ _ ZT-l JP 

.r o,_r -|- _r x_r Q — n Q r 

that is absent in a classical description. 

We remark that the projected position operators ()2.4j) differ from the Newton- Wigner 
position operators introduced in |NW49j . Whereas the latter also respect the splitting of 
the Hilbert space into particle and anti-particle subspaces, they arise from the standard 
position operator by unitary transformations. The Newton- Wigner operators are unique 



in the sense that they possess certain natural localisation properties Wig62 . However, it 
is well known that in relativistic quantum mechanics no position operators exist that leave 
the particle and anti-particle subspaces invariant, share the natural localisation properties 
and do not violate Einstein causality, see e.g. the discussion in [Tha9 2 3 . Thus even for 
free particles no complete quantum-classical correspondence exists. The goal we want to 
achieve in the following hence is to promote a receipt that allows to separate particles from 
anti-particles in a semiclassical fashion: The Hilbert space is split into mutually orthogonal 
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subspaces and observables are projected to these subspaces. Within these subspaces one 
can then employ semiclassical techniques and set up quantum-classical correspondences. 
In the next section we are going to extend this concept to particles interacting with ex- 
ternal fields. In this situation even a separation of Jri? into subspaces that are associated 
with particles and anti-particles, respectively, can only be achieved asymptotically (in the 
semiclassical limit). 



3 Semiclassical projection operators 

We now consider the Dirac equation 

d 

ih—i)(t,x) = Hif>(t,x) 

for a particle of mass m and charge e coupled to external time-independent electromagnetic 
fields E(x) = — grad0(cc) and B(x) = rot-4.(cc). The Hamiltonian therefore reads 

'h_ e 



H = ca- (-V - ~^0*0) + fimc 2 + etf>{x) . (3.1) 



For convenience we restrict our attention to smooth potentials, in which case the operator 
()3.1|) is known to be essentially self-adjoint on the domain C^°(1R 3 ) (g> C 4 , see e.g. |Tha92 j , 
and hence defines a unitary time evolution 



U(t) := e _ t 

In order to separate particles and anti-particles also in the presence of electromagnetic 
fields one would like to construct projection operators analogous to ()2.3j) that commute with 
the Hamiltonian. This, however, is not possible as can be seen by considering operators 
in a phase-space representation. Here we choose the Wigner-Weyl calculus in which an 
operator B on J^ 5 is defined in terms of a function B(x,p) on phase space taking values 
in the 4x4 matrices, 

(B^)(x) = -i— / / e ivi*-v) B (^±y^\ ^ {y) dydp . (3. 2) 

This operator is well defined on Dirac-spinors V> G y(M. 3 ) ® C 4 , if the matrix valued 
Weyl symbol B(x,p) is smooth. For this and further details of the Weyl calculus see 
|Fol89L lRoE87l IDS99j . Introducing then, e.g., the symbol 

H(x,p) = ca • (p- ~ c A { x )) + P mc2 + e <f>( x ) , ( 3 - 3 ) 

one can represent the Hamiltonian (j3.1|) as a Weyl operator (|3.2j) . For each point (x,p) in 
phase space the symbol ()3.3j) is a hermitian 4x4 matrix with the two doubly degenerate 

eigenvalues 

h±(x,p) = e<f)(x) ± y/(cp- eA(x)) 2 + m 2 c 4 . (3.4) 
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These functions can immediately be identified as the classical Hamiltonians of particles 
and anti-particles, respectively, without spin. Associated with these eigenvalues are the 
projection matrices 

2. V v(cp - eA(x)) 2 + m 2 c 4 J 

onto the respective eigenspaces in C 4 . In the absence of potentials the matrix valued 
functions ()3.5j) are independent of x and their Weyl quantisations according to ()3.2j) yield 
the projectors (|2.3|) . One could hence be tempted to view the Weyl quantisations ITq of 
the symbols (|3.5|) as appropriate substitutes for the operators ()2.3|) in the general case. 
Indeed, if the partial derivatives of A(x) are bounded, the quantities ()3.5|) and all their 
partial derivatives are bounded smooth functions of (x, p) such that the operators IIq are 
self-adjoint and bounded on the Hilbert space J$? |CV71j . The fact that the symbols (j3.5j) 
depend on x, however, implies that rig are not projection operators, but satisfy [EW96J 

(n±) a -n± = o(ft). (3.6) 

Moreover, IIq do not commute with the Hamiltonian. 

The error on the right-hand side of (|3.fij) can be improved by adding a suitable term 
of order H to the symbol 1)3. 5 j) . To guarantee that the Weyl quantisation of this symbol 
leads to a bounded operator, we now demand that the potential <fi{x) and all its partial 
derivatives are bounded, and that A(x) grows at most like some power \x\ . In this case 
the quantity 

e(x,p) := a/ (cp — eA(x)) 2 + m 2 c 4 (3.7) 

may serve as a so-called order function for the symbol ()3.3p and the construction of [BG04J 
applies. After quantisation of the corrected symbol this leads to an operator that is a 
projector up to an error of order h 2 . This procedure can be repeated to an arbitrary order 
h , see [EW96, B N99] . One hence obtains almost projection operators Il ± that almost 
commute with the Hamiltonian, 

(ft ± ) 2 - ft* = 0(h°°) and [H, ft*] = 0(/f°) , (3.8) 

meaning that the operator norm of the above expressions is smaller than any power of h. 
Furthermore, 

n+fr = o{h°°) and n + + rr = \ x + o{h°°) . (3.9) 

As a consequence of ft 1 * 1 being almost projectors, their spectrum is concentrated around 
zero and one. The standard Riesz projection formula therefore allows to construct the 
genuine orthogonal projectors 



— / (n ± -A)~ 1 dA 
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which also almost commute with H and fulfill (|3.9|) . 

Since the above construction is based on the separation into particles and anti-particles 
on a classical level ()3.4jl and is then extended order by order in h, one can view the sub- 
spaces Jif ± := P ± Jif as being semiclassically associated with particles and anti-particles, 
respectively. Moreover, due to the relation (j«18j) these subspaces are almost invariant with 
respect to the time evolution generated by H, i.e. 

U(t) P ± ^)-P ± U(t) ip = 0(th°°) for all i/j G 3? . 

Hence, up to a small error every spinor in one of the subspaces ffl^ remains under the 
time evolution within this subspace; and this is true for semiclassically long times t <C ft, , 
where N can be chosen arbitrarily large. This result is in agreement with the (heuristic) 
physical picture that particles interact with anti-particles via tunneling; the latter is a 
genuine quantum process with pair production and annihilation rates that are exponentially 
small in h. Related to this observation is the fact that eigenspinors of the Hamiltonian can 
(only) almost be associated with particles or anti-particles: If ip n G is an eigenspinor, 
Hip n = E n ip n , its projections P ± ip n are in general only almost eigenspinors (quasimodes), 
i.e. 

(H-E n )P ± Tp n = 0(h°°) . 

Thus the discrete spectrum of the Hamiltonian cannot truly be divided into a particle and 
an anti-particle part. Further consequences will be investigated in section |3 

The above discussion raises the question whether the semiclassical projectors P ± can 
(at least semiclassically) be related to spectral projections of the Hamiltonian. Accord- 
ing to (|2.3|) . in the case without potentials the operators Pq are obviously the spectral 
projectors to (— oo, — mc 2 ) and (mc 2 , oo), which are exactly the (absolutely continuous) 
spectral stretches associated with particles and anti-particles, respectively. In the presence 
of potentials, depending on their behaviour at infinity, there exist constants E + > £L such 
that the spectrum inside (£L, E + ) is discrete and absolutely continuous outside; e.g., if the 
potentials and their derivatives vanish at infinity one finds E± = ±mc 2 , see |ThaQ2| . For 
E G (-&_, E + ) not in the spectrum of H we then denote the spectral projectors to (—00, E) 
and (E, 00) by Pg and P£, respectively. For the semiclassical considerations to follow we 
also assume that 

sup h-(x,p) < E < inf h + (x,p) ; 

in particular, it is required here that the classical particle and anti-particle Hamiltonians 
are separated by a gap as (x, p) runs over phase space. Therefore, the two classical energy 
shells 

M% = {(x,p); h±(x,p) = E} 

must not intersect nor touch. On the submanifold p = -A(x) this condition requires that 
the variation of the potential 4>{x) is strictly restricted by 2mc 2 . In such a situation the 
semiclassical projectors P^ 1 are semiclassically close to the spectral projectors, 

||P± _P±|| = 0(fr°°) . 
A full proof of this statement can be found in [BG04 . 
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4 Semiclassical observables 



The example of the standard position operator for the free particle and its Zitterbewegung 
demonstrates that not all quantum observables possess a direct (semi-) classical interpre- 
tation. Only the diagonal blocks with respect to the projections ig allow for establishing 
quantum-classical correspondences. This observation applies in particular to dynamical 
questions because only after projection there exists an unambiguous classical Hamiltonian 
(pOJ. Following the discussion in the previous section, these remarks carry over to the 
general case appropriately. Our aim therefore is to identify a sufficiently large class of 
observables with a semiclassical meaning. We moreover require this class to be invariant 
with respect to the quantum time evolution, i.e., B being of this class should imply that 
B{t) = U(t)*BU(t) is of the same type. 

A natural starting point in the present context is to choose the Weyl representation 
(IH.2j) for operators on M' . It is well known that this is possible for a fairly large class of 
operators, namely those that are continuous from ^(M 3 x M. 3 ) Cg>C 4 to ^'(IR 3 x M. 3 ) ®C 4 , see 
e.g. |Fol89tlDS99j . In order to achieve manageable algebraic properties, however, one must 
restrict the class of operators further. For convenience we here choose Weyl quantisations 
of smooth symbols B(x,p) that are, along with all of their partial derivatives, bounded 
on K 3 x Kl This restriction leads to bounded operators on J^, see |CV71| . We admit 
^-dependent symbols, since even when restricting to quantisations B of ^-independent 
symbols B(x,p) their time evolutions B(t) will necessarily acquire an ^-dependence in 
their symbols B(t)(x,p; h). But for the purpose of convenient semiclassical asymptotics 
we restrict the ^-dependence in a particular way, 

oo 

B(x,p;h)~J2fr kB k&P) • (4- 1 ) 

fc=0 

The asymptotic expansion is to be understood in the sense that 

jV-1 

\%a;(B{x,nh)-^h h B k (x,pj)\\ 4 <c^h N (4.2) 

fc=0 

holds for all (x,p) G M 3 x M 3 and all iV G N uniformly in h as h — > 0, where || • \\4x4 
denotes an arbitrary matrix norm. We call the quantisations of such symbols semiclassical 
observables. These form a subalgebra O sc in the algebra of bounded operators on M' . As 
usual, only the self-adjoint elements of O sc are true quantum mechanical observables; they 
can be characterised by symbols taking values in the hermitian 4x4 matrices. 

The crucial point now is to identify a subalgebra O im in O sc that is invariant with 
respect to the quantum time evolution. We recall that in the case of scalar operators, such 
as Schrodinger-Hamiltonians on L 2 (M 3 ), the respective algebra O sc itself is invariant, see 
e.g. [Rob87j. The Zitterbewegung of a free relativistic particle, however, shows that in the 
present context this can no longer hold. This example also suggests that O uw might consist 
of those semiclassical observables that are block-diagonal with respect to the semiclassical 
projections 11^. In the following we are going to demonstrate that this is indeed the case. 
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Since the necessary constructions take place in terms of symbols, we first recall how 
the algebraic properties of semiclassical observables are reflected on this level, see e.g. 
[Fol89, Rob87, DS99 : The product of two semiclassical observables B,C G O sc is a Weyl 
operator with a symbol denoted B#C(x,p; h). This symbol has an asymptotic expansion 
of the type (jO|l - (fOjl that reads 



j,k,l=0 



y=x 



We remark that if the operator C is replaced by the Hamiltonian H a corresponding product 
formula for the symbol B#H(x,p; H) exists, with the exception that on the right-hand side 
of the estimates ()4.2j) the order function ()3.7|) appears as an additional factor. 

We are now in a position to write down the Heisenberg equation of motion for the time 
evolution B(t) of a semiclassical observable in terms of symbols, 

^B(t)(x,p;h) = ^_(H#B(t) - B(t)#H)(x,p;h) . (4.3) 

The strategy we follow now is to assume an asymptotic expansion of the form (|4.2j) for 
the symbol B(t)(x,p; H), insert this into (|4.3JI and solve, if possible, the hierarchy of equa- 
tions for the coefficients B(t)k(x,p) of h k . A first obstacle to this procedure arises upon 
investigating the leading terms on the right-hand side of f)4.3j) . 

^B(t) = l -[H,B{t) \ - ±({H,B{t) ] - {Bit^H^+iiH^m+Oih) . (4.4) 

Here [•, •] is a matrix commutator and 

{B,C}(x,p) = (V P B- V X C- V x B-V b C)(x,p) 

denotes the Poisson bracket for matrix valued functions on phase space. Due to the factor 
of 1/h in the first term on the right-hand side of (|4.4|) the proposed strategy can only 
be consistent if the matrix valued function B(t) (x,p) is block- diagonal with respect to 
the eigenprojections n o (x,p) of the symbol matrix H(x,p). This condition is the lowest 
semiclassical order of the expected restriction on the operators in O mv . 

Necessary and sufficient conditions to be imposed on semiclassical observables to be in 
O irw can be obtained in any semiclassical order, if the solvability of the equations (|4.3jl for 
the diagonal and for the off-diagonal blocks of B(t) with respect to the projection symbols 
IP 1 are investigated. This leads to: 

Theorem 4.1. A semiclassical observable B G O sc lies in the invariant algebra O mv , if and 
only if its off-diagonal blocks with respect to the semiclassical projections P are smaller 
than any power of h, 

P ± BP T = 0(h°°) . (4.5) 
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This statement is a version appropriate for the Dirac equation of a result contained in 
[BG04j. There one can also find a proof. 

The Hamiltonian (j3.1j) is not in O imr) but this is only due to our restriction to bounded 
operators; otherwise, the relations (pEE} and imply that H indeed possesses the prop- 
erty ()4.5j) . Furthermore, up to errors of the order h°° the time evolution of the diagonal 
blocks P ± BP ± is governed by the projected Hamiltonians /f P . An interpretation of The- 
orem now is rather obvious: An observable remains semiclassical under the quantum 
dynamics if it does not contain off-diagonal blocks representing an interference of particle 
and anti-particle dynamics, since this has no classical equivalent. Terms exceptional to 
this must be smaller than any power of h, indicating that they arise from pure quantum 
effects. 

5 Classical limit of quantum dynamics 

Given an observable B e O mY , we are interested in the connection between its time evolu- 
tion B(t) and appropriate classical dynamics. The diagonal blocks of B are (approximately) 
propagated with the projected Hamiltonians HP ± , whose symbols read 

H(x,p) I${x,p) = h ± (x,p) I^(x,p) 

to leading order. We therefore expect the equations of motion generated by the eigenvalue 
functions (jH.4|) . 

x±(t) = {h ± ,x ± (t)} , p ± {t) = {h ± ,p ± {t)} , (5.1) 

to play an important role for the classical limit of the quantum time evolution. The solu- 
tions to ()5.1|) define the Hamiltonian flows $±(p, x) = (p±(t), x±(t)), with (p ± (0), x±(0)) = 
(p, x), on phase space. These flows represent the classical dynamics of relativistic particles 
(and anti-particles). Spin is absent from these expressions, reflecting the fact that a priori 
spin is a quantum mechanical concept and its classical counterpart has to be recovered in 
systematic semiclassical approximations of the quantum system. 

In the present situation the kinematics and the dynamics of spin are encoded in the 
matrix part of the symbols of observables. For an observable B e O imr the leading semi- 
classical order of its symbol is composed of the functions 

(l$B I$)(x,p) (5.2) 

taking values in the (hermitian) 4x4 matrices. These matrices act on the two dimensional 
subspaces n^(:r,p)C 4 of C 4 , which can be viewed as Hilbert spaces of a spin 1/2 attached 
to a classical particle or anti-particle, respectively, at the point (x,p) in phase space. In 
order to work in these two dimensional spaces explicitly, one can introduce orthonormal 
bases {e^(x,p), ef(x,p)} for each of them. For convenience the vectors e^(x,p) E C 4 
are chosen as eigenvectors of the symbol matrix H(x,p) with eigenvalues h±(x,p). Upon 
expanding vectors from n^(:r,p)C 4 in these bases one introduces isometries V±(x,p) : 
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C 2 -> rhf (a;,»)C 4 , such that V±(x,p)V£(x,p) = I^(x,p) and V£(a:,p)V ± (a:,p) = 1 2 - A 
possible choice for these isometries is 

K f a: = 1 ( fcfo P) + mc ') M 

+ v /2£(;r,p)(£(a ; ,p)+mc 2 ) \ (cp - eA(x)) - cry ' 

^ / v = 1 ( (cp - eA(x)) ■ <x \ 

" v /2 £ ( ; r,p)(£(a ; ,p)+mc 2 ) V~ (efo p) + ™c 2 ) t 2 ) 



That way the diagonal blocks ()5.2j) can be represented in terms of the 2x2 matrices 

(VlB V±)(x,p) , (5.3) 

which are hermitian if Bq(x,p) is hermitian on C 4 . 

The matrix valued functions ()5.3|) can be mapped in a one-to-one manner to real valued 
functions on the sphere S 2 via 

b£(x,p,n) := ^tr(A 1/2 (n) (V£B Q V ± ){x,pj) , (5.4) 

where n E R 3 with \n\ = 1 is viewed as a point on S 2 and 

A 1/2 (n) := i(l 2 + V^n-cr) 

takes values in the hermitian 2x2 matrices, see (VGB89j. The quantisers Ai/ 2 (n) provide 
a quantum-classical correspondence on the sphere that is covariant with respect to SU(2)- 
rotations in the following sense, 

g A 1/2 (n) g- 1 = A 1/2 (R(g)n) . (5.5) 

Here R(g) E SO(3) is a rotation associated with g E SU(2) through gn-ag" 1 = (R(g)n)-a, 
see |BGK01j for further details. The semiclassically leading term of a (block-diagonal) 
semiclassical observable B can therefore be represented by two real valued functions on 
the space R 3 x K 3 x S 2 . Through the relation 

{V±B V±)(x,p)= I b$(x,p,n)A 1/2 (n) dn , 
Js 2 

where dn is the normalised area form on S 2 , the matrix valued expressions ()5.3|) are un- 
ambiguously recovered from the functions (|5.4|) [BGK01 . 

The point n E S 2 that arises as an additional variable in a general classical observable 
(15 .4|) can be viewed as a classical equivalent of spin. This interpretation is suggested by 
the fact that the classical observable ()5.4|) associated with quantum spin, 

B — — Ei. with Et. := \ ^ ^ 
2 \0 a k 
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can be calculated as 

bo (x,p,n) = ^ tr^A 1/2 (n) V±{x,p) Z k V±(x,pfj 

= v / l(l + 1 )^^ • 

Up to its normalisation, which depends on h and the spin quantum number s = 1/2, the 
point n G S 2 can thus be considered as a classical spin. The latter is therefore represented 
on its natural phase space S 2 , and the space I 3 x K 3 x S 2 on which the classical observables 
f!5.4|) are defined can be viewed as the combined phase space of the translational and the 
spin degrees of freedoms. 

We now want to identify the classical dynamics of both the translational and the spin 
degrees of freedom that emerge in the semiclassical limit of the quantum dynamics gen- 
erated by the Hamiltonian H. To this end we recall that 0\ m was defined to contain 
those semiclassical observables B G O sc whose time evolutions B(t) = U(t)*BU(t) are 
again semiclassical observables. Hence B(t) is a Weyl operator of the form ()3.2|) with sym- 
bol B(t)(x,p; h) that possesses an asymptotic expansion of the type (|4.1|) . The leading 
term B(t)o(x,p) in this expansion then yields the classical time evolution of both types of 
degrees of freedom. 

Theorem 5.1. When B G 0- mv the semiclassically leading term in the asymptotic expan- 
sion of the symbol of its time evolution B(t) reads 

B(t) (x,p)= dl(x,p,t)(I%B I%)(& v (x,p))d u (x,p,t) . (5.6) 

^e{+,-} 

The unitary Ax A matrices d±(x,p,t) are determined by the transport equations 

d ± (x,p,t)+iH ± (&±{x,p))d ± (p,x,t) = , d ± (p, x, 0) = 1 4 , (5.7) 
where the effective spin-Hamiltonians H± are defined as 

h ± : = -i[n± {h±, n±>] + ~ (/> ± n±{n±, n±}n± + n ±{n± h - /i ± n ±}n±) . (5.8) 

We refrain from giving explicit expressions for the 4x4 effective spin-Hamiltonians ()5.8|) 
here since below we will only work in a 2 x 2 representation; instead we refer to |Spo00| . A 
statement of the type made in Theorem 15.11 is usually called an Egorov theorem [Ego69 . 
The present version is covered by the Egorov theorem in IKK) lj. where one can also find 
a proof. 

The 4x4 matrices d±(x,p,t) can be shown to map the two dimensional subspaces 
n^(cc,p)C 4 of C 4 unitarily to the propagated subspaces n^($' ± (a;, p))C 4 |BG04| . They 
hence describe the transport of the (anti-) particle spin along the classical trajectories 
(x±(t) : p±(t)). If one prefers to work in the orthonormal eigenbases {ef (x, p), e^(x, p)} 
of the spaces n^(a;,p)C 4 one can introduce the unitary 2x2 matrices 

D ± (x, p, t) := VI (¥ ± (x, p)) d ± (x, p, t) V±(x, p) , 
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which are determined by the equations 

D ± (x,p,t)+ 1 -C ± ($i(x,p)) ■crD ± (p,x,t) = , D ± (p, x, 0) = 1 2 , (5.9) 

following from ()5.7|) . The transformed effective spin-Hamiltonians C± ■ er/2 derive from 
()5.8|) and can be expressed in terms of the electromagnetic fields through 

where e(x,p) is defined in f)3 . 7j) . Since therefore the effective spin-Hamiltonians are her- 
mitian and traceless 2x2 matrices, the solutions D± of (|5.9|) are in SU(2). Related results 
have previously been obtained in a WKB-type situation [IRK631 IBK981 lBK9_9J and in the 
context of a semiclassical propagation of Wigner functions |Spo00 . 

The transport matrices D±(x,p,t) carry two-component spinors along the trajectories 
$j_(£C,p) and induce a classical spin dynamics along (anti-) particle trajectories. These 
combined classical dynamics can be recovered upon representing the leading symbol (|5.6J) 
in terms of the functions (|5.4jl on the combined phase space, 

b(t)t(x,p,n)=tr(A 1/2 (n) (V* B(t) V ± ){x,p)) 

= tv(D ± (x,p,t) Ai /a (n) D* ± (x, Pl t) (V± B V ± ) ($* ± (aj,p))) 
= #($' ± (a:,p),i2p ± (:c,p,t))n) , 

where in the last line the covariance ()5.5|) has been used. Hence, in leading semiclassical 
order the dynamics of an observable B e O inv can be expressed in terms of the classical 
time evolutions 

(x,p) h-> <&±(a:,p) and n i-> n±(i) = R(D±(x,p,t))n . (5.10) 

The spin motion on the sphere thus emerging obeys the equation 

n±(t) = C±($i(x,p)) x n±(t) , n±(0) = n , (5.11) 

that is implied by (|5.9j) . These spin dynamics exactly coincide with the Thomas precession 
that was derived in a purely classical context |Tho27j . 

Both types of time evolutions in (|5.10J) can be combined to yield the dynamics 

(x,p,n)^Yl(x,p,n) := (& ± (x,p),R(D ± (x,p,t))n) (5.12) 

on the phase space R 3 x I 3 x S 2 . The properties 

Yl=°{x,p,n) = (x,p,n) and Y* +t (x,p,n) = Y*(Yl(x,p,n)) 

can easily be established, showing that ()5.12|) defines flows on the combined phase space. 
These flows are composed of the Hamiltonian flows defined in (|5.1|) on the ordinary 
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phase space M 3 x R 3 and the spin dynamics (|5.11|) on S 2 that are driven by the Hamiltonian 
flows. Dynamical systems of this type are known as skew-product flows, see e.g. [CFS82 . 
They leave the normalised measures 

d£%(x,p) dn = — ^—c5(h±(x,p) — E) dxdpdn 
volll E 

on the combined phase space invariant that are products of Liouville measure d£^(x,p) 
on (microcanonical ensemble) and the normalised area measure dn on S 2 . We hence 
now conclude that the classical limit of the quantum dynamics generated by the Dirac- 
Hamiltonian (j3.1j) is given by the two skew-product flows (j5.12|) combining the Hamiltonian 
relativistic motion of (anti-) particles with the spin precession along the (anti-) particle 
trajectories. 



6 Semiclassical behaviour of eigenspinors 

In section E] we saw that in general eigenspinors ip n of the Dirac-Hamiltonian (|3.1|) cannot 
uniquely be associated with either particles or anti-particles. For the purpose of semi- 
classical studies we therefore prefer to work with the normalised projected eigenspinors 



IIP^II (6 ' 1) 

Due to the fact that the projectors almost commute with the Hamiltonian (|3.8|) . the quan- 
tities (jfi.lj) are almost eigenspinors (quasimodes) of both H and HP ± , 



\\(H-E n )<t>5\\=r5 and \\(HP ± -E n )^\\ =s± , 
with error terms and given by 

± WlH^P^nW , ± \\[hp ± ,P ± }A, 

r„ = — — - and sf " 



Thus, if the norms \\P ip n \\ of the projected eigenspinors are not too small, i.e. if ||P ip n \\ > 
Ch N holds for some C > and N < oo, the error terms and are of size h°°. 

This property implies that if the spectra of the operators H and HP ± are discrete 
in the intervals [E n — rj, E n + r*] and [E n — , E n + s^], respectively, these operators 
each possess at least one eigenvalue in the respective interval, see e.g. |Laz93j . For H the 
statement is trivial since E n is known to be an eigenvalue, but regarding HP ± it provides 
new insight. In particular, the mean spectral densities of the operators HP ± are in general 
smaller than that of H, see (16 .2|) below. Hence, the eigenvalues E n of H must cluster 
in such a way that successive intervals [E n — r^,E n + r^] overlap and therefore several 
such intervals share eigenvalues of HP ± . That way the separation of eigenvalues of H 
is linked with the norms of projected eigenspinors P ± ip n . Suppose, e.g., all norms were 
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bounded from below by H-F^^H > CH N , such that = 0(h°°), then sufficiently many 

eigenvalues E n have to cluster on a scale 0(h°°). On the other hand, if the separations 
of neighbouring eigenvalues of H were known to be bounded from below by h , a certain 
number of eigenspinor-projections would have to possess norms of size h°° in order to 
produce sufficiently large errors r*. 

More quantitative statements can be made on the ground of Weyl's law for the numbers 
N EjUJ and N Euj of eigenvalues the operators H and HP ± , respectively, possess in the interval 
[E — hu), E + flu)]. Namely, as h —>■ 

2uj vol Q E + vol Q E 

~ " ^ ' (6.2) 
i_ 2tu volf2p 



Thus, if at energy E both energy shells fi^ in phase space have positive volumes, the ratios 
u)/Ne,u governing the previous discussion are semiclassically determined by the relative 
fraction of the volumes of the associated energy shells. Moreover, a version of Weyl's law 
that includes expectation values of operators (Szego limit formula) yields [BG04 

N%„~ \\ p± ^\\ 2 - (6.3) 

E-hw<E n <E+riuj 

This relation may suggest two (extreme) scenarios: 

(i) Roughly iVj 1 ^ of the projected eigenspinors P ± ip n are close to ip n and the remaining 
projections are semiclassically small. 

(ii) The projections P ± ip n equidistribute on both energy shells in the sense that each 

VVi|| 2 is approximately w /Ne,u- 

On the ground of the existing knowledge one cannot exclude either alternative, nor a 
mixture of both. However, a comparison with symmetric vs. asymmetric double-well 
potentials, where an eigenfunction either localises in one well (asymmetric case, compare 
(i)), or has noticeable projections to both wells (symmetric case, compare (ii)), suggests 
that in generic cases scenario (i) should be expected. 

What is possible, though, is to estimate the number of projected eigenspinors with 
norms that are bounded from below, 

Nf := #{n; E - hw < E n < E + hw, HP^nf > $} ■ 

From (JH2J) and (JEHJ) it follows that 

lim > t - 5 . 

h-+o Ne,w vol Vl E + vol Vl E 
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Therefore, if S is sufficiently small (independent of K) the right-hand side is positive and 
thus a finite portion of the projected eigenspinors have semiclassically non-vanishing norms. 
This observation allows to take the associated normalised projected eigenspinors (JHU into 
account for the following considerations. 

We now explore the consequences an ergodic behaviour of the classical skew-product dy- 
namics Y± exerts on the projected eigenspinors. In an analogous situation for Schrodinger- 
Hamiltonians the ergodicity of the associated Hamiltonian flow implies that the Wigner 
transforms of almost all eigenfunctions equidistribute on the corresponding energy shell. 
This result, known as quantum ergodicity, goes back to Shnirelman [Shn74j and has been 
fully proven in |Zel87[ ICdV85[ IHMR87J . In the case under study we now suppose that the 
energy E lies in an interval (E_,E + ), in which the spectrum of the Dirac-Hamiltonian is 
discrete. Moreover, at least one of the classical energy shells and shall be non-empty 
and the periodic orbits of the Hamiltonian flows $± shall be of measure zero on f2^. If 
then Y± is ergodic on Q E x S 2 , the time average of a classical observable of the type (|5.4jl 
equals its phase-space average, 

lim 1 [ T b±(Yi(x,p,n)) dt= I [ b±(x',p',n') d£%{x\p') dri =: M±(fo±) , 
T ^°° J- Jo Ja% Js 2 

for almost all initial conditions (x,p,n). In this setting the main result on quantum 
ergodicity of projected eigenspinors is (for a full proof see [B G04j ): 

Theorem 6.1. Under the conditions stated above, in particular if the skew-product flow Y± 
is ergodic on fi^ x S 2 , in every sequence {<^} n eN of normalised projected eigenspinors with 
associated eigenvalues E n 6 [E — Jtlu, E + tkj] and \\P ± ip n \\ 2 > 5 there exists a subsequence 
{(ptjam of density one, i.e. 

lim #{q; \\P^ n J 2 >5} l 
#{r2; ||P±Vn|| 2 >5} 

such that for every semiclassical observable B e O sc 

hm(^ Q , J B0^) = M|(6 o ± ) . (6.4) 

The density-one subsequence {0„ } a eN can be chosen independent of the observable. 

The statement of this theorem, that in case of classical ergodicity quantum mechanical 
expectation values converge to a classical ergodic mean, can be rephrased in terms of 
more explicit phase-space lifts. To this end one first introduces the matrix valued Wigner 
functions 

W[ip)(x,p) := / e~Ti pv ^j{x- \y) ®^{x + \y) Ay 

associated with any ip 6 L 2 (M 3 ) ® C 4 . Following the prescription (|5.4|) this 4x4 matrix 
can be converted into scalar form through 

w±[i>](x,p,n) := itr(A 1/2 (n) (V£ Wty] V±)(x,pj) . 
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Due to the projections inherent in the expectation value on left-hand side of (J6.4J) only 
one diagonal block of the observable contributes, so that this expectation value can be 
reformulated as 



' II:, ' ^ 



j^- 3 J j tr(W[0± }P ± #B#P ± )(x 1 p) dxdp 
—J— JJ tT^W[^ a ]V ± )(V^B V ± )(x,p)(l + O(h))) dxdp 
{w±[0nJ( x 'Pi n ) b o( x iPi n )( 1 + °( h ))) dxdpdn . 



1 




(2nh) 3 

The principal result ()6.4|) of quantum ergodicity can hence be read as saying that 
lim ," U ' = T 5(h±(x,p)-E) . 

This relation has to be understood in a weak sense, i.e. after integration with a symbol 
b^(x,p,n). Hence, the scalar Wigner transforms of projected eigenspinors become semi- 
classically equidistributed on the associated phase space x S 2 once the classical time 
evolution Y± is ergodic on this space. 

As the discussion at the beginning of this section shows, the difficulties with statements 
about genuine eigenspinors derive from the coexistence of the particle and anti-particle 
subspaces; and because interactions between these subspaces on a scale h°° cannot be 
controlled within the present setting. However, if one of the energy shells were empty 
there is only one classical manifold onto which phase-space lifts of eigenspinors could 
condense, namely x S 2 . In such a case, say when is empty, the statement of 
Theorem 16.11 applies to a density-one subsequence {^ n J of eigenspinors themselves, such 
that 

lire y (ij na ,P+BP+ij n J = M±(fo±) and lim p-Bp-^J = . (6.5) 

A corresponding statement then holds for the associated scalar Wigner transforms w±[i/j na }. 
It has been mentioned previously that the restriction of having only one non-empty energy 
shell at energy E requires potentials that do not vary too strongly; e.g., <j>(x) must not 
vary as much as 2mc 2 . This condition is fulfilled in many physically relevant situations, in 
which (|6.5J) therefore applies if only the classical particle dynamics Yf are ergodic. 
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